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We study numerically the geometrical properties of minimally weighted paths that appear in the 
negative-weight percolation (NWP) model on two-dimensional lattices assuming a combination of 
periodic and free boundary conditions (BCs). Each realization of the disorder consists of a random 
fraction (1— p) of bonds with unit strength and a fraction p of bond strengths drawn from a Gaussian 
distribution with zero mean and unit width. For each such sample, the path is forced to span the 
lattice along the direction with the free BCs. The path and a set of negatively weighted loops 
form a ground state (GS). A ground state on such a lattice can be determined performing a non- 
trivial transformation of the original graph and applying sophisticated matching algorithms. Here 
we examine whether the geometrical properties of the paths are in accordance with predictions of 
Schramm-Loewner evolution (SLE). Measuring the fractal dimension and reviewing Schramm's left 
passage formula indicates that the paths cannot be described in terms of SLE. 
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I. INTRODUCTION 

The statistical properties of lattice-path models on 
graphs, often equipped with quenched disorder, have ex- 
perienced much attention during the last decades. Such 
paths can be as simple as boundaries of percolation clus- 
ters p] or domain walls in planar Ising systems [2]. Be- 
yond this, they have proven to be useful in order to de- 
scribe more complex line-like quantities as, e.g., linear 
polymers in disordered media [3HS]) vortices in high T c 
superconductivity pj [8] , cosmic strings in the early uni- 
verse [SHU]; and domain- wall excitations in disordered 
systems such as 2d spin glasses [T2J H3] and the 2d solid- 
on-solid model [L4j . So as to analyze the statistical prop- 
erties of these lattice path models, geometric observables 
and scaling concepts similar to those used in percola- 
tion theory [TJ [T5] or other string-bearing models [151 HZ] 
are often applicable. Recently, the question which of 
these path models can be described in terms of Schramm- 
Loewner evolution (SLE) [T81E0] , has attracted much at- 
tention. 

Basic problems like percolation or domain walls in 
Ising systems can be treated numerically by standard ap- 
proaches [TJ [2T]. For more complex models, this is not 
straight-forward. Nevertheless, the precise computation 
of these paths can often be formulated in terms of a com- 
binatorial optimization problem and hence might allow 
for the application of exact optimization algorithms de- 
veloped in computer science [22] [23] . 

In this article we study paths in the negative-weight 
percolation (NWP) model [23]. The basis of this model 
is given by a weighted, undirected, d-dimensional lattice 
graph with side length L where weights are assigned to 



'Electronic address: christoph.norrenbrock@uni-oldenburg.de 
^Electronic address: olivcr.mclchert@uni-oldenburg.de 
tElectronic address: alexander.hartmann@uni-oldcnburg.de 



L 




FIG. 1: Illustration of minimum-weight configurations con- 
sisting of loops (gray) and one path (black) in a regular 2d 
square lattice of side length L = 64. The boundary conditions 
are periodic in horizontal (left-right) direction and free in ver- 
tical direction. The path is forced to connect the middle of 
the lower boundary to the upper boundary. The minimum- 
weight configurations minimize Eq. j3J. For a small value of 
the disorder parameter p = 0.2, the path crosses the lattice 
in a rather straight fashion, since the path weight is mainly 
determined by its length. Increasing the order parameter to 
p = 0.3413, the path length increases as well and the path 
takes advantage of more negative-weighted edges. If the dis- 
order parameter is increased to p = 0.6, the loops and the path 
will become more dense. 



the edges taken from a distribution that allows for values 
of either sign. For a given realization of the disorder, a 
configuration consisting of one path, spanning the lattice 
along the direction of the free boundaries (the path ex- 
tremities are forced to lie on the free boundaries), and 
possibly a set of loops, i.e. closed paths, can be found 
such that the loops and the path do not intersect one 
another and the total sum of the weights assigned to the 
edges which build up this configuration attains an ex- 
act minimum. The fraction of negatively weighted edges 
can be controlled by a tunable disorder parameter p that 
gives rise to qualitatively different ground states as de- 
picted in Fig. [I] As presented here, the NWP problem is 
a theoretical model of intrinsic interest. As an example 
that makes use of the statement of the NWP problem, 
an agent can be imagined that travels on a graph. While 
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traversing an edge, the agent either needs to pay some 
resource (signified by a positive edge weight) or is able 
to harvest some resource (signified by a negative edge 
weight). In order to gain as many resources as possi- 
ble, the path/loops obtained in the context of the NWP 
problem serve as a guide to find optimal routes. 

The problem of detecting the exact minimum-weight 
configuration for a given realization of the disorder be- 
comes solvable through a mapping of the original graph 
to the minimum-weight perfect-matching optimization 
problem as outlined in Section [n] This mapping allows 
for applying exact polynomial-time-running algorithms, 
thus large instances can be solved. Note that the same 
mapping and algorithms can be used in the context of 
finding ground-state spin configurations for the planar 
triangular random-bond Ising model, since this problem 
is equivalent to the NWP problem on a planar honey- 
comb lattice |25ll26]. 

In a previous study it was shown that the NWP model 
features a disorder-driven, geometric phase transition 
[24] . Depending on the disorder parameter, two distinct 
phases can be identified: (i) a phase where the loops (if 
there are any) are small and the shape of the path is 
rather straight-lined, reflecting a self-affine scaling of the 
path length (cf. Fig. [IJa)) and (ii) a phase where large 
loops emerge and the extension of the horizontal projec- 
tion of the path is O(L), reflecting a self-similar scaling of 
the path length (cf. Fig.[ljc)). In the limit of large system 
sizes, there is a particular value of the disorder param- 
eter at which the path simultaneously spans the lattice 
along both directions for the first time (cf. Fig. |ljb) ) . 
For two-dimensional (2D) lattice graphs, the respective 
disorder-driven phase transition was investigated using 
finite-size scaling analyses and characterized by a set of 
critical exponents. It turned out that the exponents were 
universal in 2D and clearly distinct from those describ- 
ing other percolation phenomena. Regarding loops only, 
the influence of dilution on the critical properties of the 
2D NWP phenomenon was investigated in a subsequent 
study [27 . In a further study, the upper critical dimen- 
sion d u of the NWP model was searched for by perform- 
ing simulations on hyper-cubic lattice graphs in dimen- 
sions d — 2 through 7 and evidence was found for an 
upper critical dimension d u —6. In contrast to that, an- 
other upper critical dimension <i^ PL = 3 was determined 
for densely packed loops which appear if the disorder pa- 
rameter is far above the critical point that indicates the 
phase transition [29]. Only recently, the mean-field be- 
havior of the NWP model was investigated on a random 
graph with fixed connectivity [30] . By means of numer- 
ical simulations as well as an analytic approach (using 
the replica symmetric cavity method for a related poly- 
mer problem), the location of the phase transition and 
the values of the critical exponents were determined. 

Here we study the geometrical properties of minimum- 
weight paths connecting opposite (free) boundaries in the 
2D NWP model at the critical point. To be more pre- 
cise, we find the answer to the question whether these 



paths belong to the particular family of planar curves 
that can be described in terms of SLE [I8H20] . The 
(chordal) SLE formulation describes a self-avoiding curve 
7t ("time" t £ ) in the upper half plane H, which is 
equivalent to any other simply connected planar domain 
by conformal invariance. The curve jt grows dynamically 
from the origin to infinity with time t. A conformal map- 
ping gt(z) (z £ C corresponding to the two-dimensional 
plane) maps the complement of the set consisting of 74 
plus all points that are not reachable from infinity with- 
out crossing j t , i.e. the hull, to the upper half plane at 
each infinitesimal time step. Thus, the boundary between 
reachable and unreachable points is mapped on the real 
axis and gt(z) maps the growing tip successively to a 
real point a t . As a consequence, a t moves continuously 
on the real axis during the growing process. Then the 
time-evolution of gt(z) satisfies an ordinary differential 
equation, also referred to as Loewner equation: 

- -TT— W 

dt g t (z)~a t 

If at is proportional to an one-dimensional Brownian mo- 
tion \/HBt, which is a stochastic Markov process, the 
resulting curve will be conformal invariant |31j and the 
evolution will be referred to as SLE. Consequently an en- 
semble of curves that can be described in terms of SLE 
is characterized and classified by just one parameter, the 
so-called diffusion constant k. In the last decade, a lot 
of physical systems have been studied in the context of 
SLE. To name but a few, boundaries of vorticity clusters 
in 2D turbulence (ft = 6) [32], spin cluster boundaries in 
the Ising model at the critical temperature (k = 3) [33] , 
and critical percolation (k = 6) [34] fall into the classifi- 
cation scheme of SLE. Numerical evidence [35] indicates 
that also two-dimensional spin-glasses can be described 
in terms of SLE. In contrast to these examples, as it will 
be shown in this article, this is not the case for the NWP 
model. 

The remainder of the presented article is organized as 
follows. In Section [Ii] we introduce the model in more 
detail and outline the algorithm used to compute the 
minimum-weight configurations consisting of loops plus 
one path. In Section |TTT} we describe two measuring tech- 
niques for determining the diffusion constant and present 
the results of our numerical simulations. In Section lTVl we 
conclude with a summary. Note that an extensive sum- 
mary of this paper is available at the papercore database 



II. MODEL AND ALGORITHM 

In this article we study 2D square lattice graphs G = 
(V, E) with both rectangular and circular shape. Never- 
theless, in this section the model and the algorithm are 
described on the basis of the rectangular form only, since 
there is no conceptual distinction between either shapes. 



Some details concerning the different shapes are outlined 
in Section UTI1 

The top and bottom boundaries of the rectangular 
square lattice graphs are chosen to be free and the left 
and right ones are periodic, i.e. the considered lattices 
are cylindrical. The system size L — L y denotes the 
number of nodes in the vertical direction. As justified in 
Section |III| the number of nodes in the horizontal profile 
corresponds to L x = AL + 1. Thus the considered graphs 
have N — \V\ — i(4L + l) sites i £ V and a number of 
\E\ = (2£-l)(4L+l) undirected edges {i,j} G E that 
join adjacent sites i, j € V. We further assign a weight 
Wij to each edge contained in E. These weights repre- 
sent quenched random variables that introduce disorder 
to the lattice and are drawn from a distribution that al- 
lows for values of either sign. In the presented article the 
edge weights are taken from a "Gauss-like" distribution 
characterized by a tunable disorder parameter p: 



(2) 
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(b) 



P{uj) = (1 - p)S(oj - 1) + pexpf-w 2 ) /V2tt, 

< p< 1. 



This weight distribution explicitly allows for loops £ with 
negative weight given by u>c = Y^u j}ec ^ n addition 
one path V is forced on the lattice whose endpoints lie on 
the horizontal boundaries. One cndpoint is fixed at the 
central node of the lower boundary and the other end ter- 
minates at a site at the top boundary determined by the 
optimality criterion given below. Moreover, the loops and 
the path are not allowed to intersect each other, hence 
every edge is contained in at most one loop (or path). 
Given a realization of the disorder, the exact geometri- 
cal configuration C of loops and one path is determined 
by the following criterion: The configurational energy £ 
defined as 



e = J2^ 

cec 



- UJ-p 



(3) 



has to be minimized, where uj-p denotes the total of all 
edge weights belonging to the path. Such minimum- 
weight configurations for different values of the disorder 
parameter are presented in Fig. [T] 

To get over that "minimum- weight configuration prob- 
lem" we transform the original graph as detailed in |37j 
to an appropriate auxiliary graph and detect a minimum- 
weight perfect matching (MWPM) gSJ [351 I3S] subse- 
quently. Here, we give a concise description of the algo- 
rithmic procedure illustrated in Fig. [2] for a given realiza- 
tion of the disorder (for clarity without periodic bound- 
ary conditions): 

(1) In order to force a path that satisfies the specified 
boundary conditions, two "extra" nodes are added in the 
following way: One extra node is linked to the node lo- 
cated in the center of the lower boundary and the other 
extra node is connected to all nodes at the upper bound- 
ary (c.f. Fig.[2jb)). All added edges get zero weight. 

(2) Each edge is replaced by two "additional" nodes 
and three edges which are arranged in a row as depicted 
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FIG. 2: Illustration of the algorithmic procedure: Illus- 
trated edges that are not marked by a weight visually, carry 
a zero weight actually (except in (d) and (e)). (a) The ini- 
tial, weighted lattice graph, (b) Two extra nodes (gray) are 
added to render a path possible. They are connected to the 
lattice boundaries as described in the text, (c) Each edge 
is replaced by two nodes (black boxes) and three edges join- 
ing the original nodes and the added nodes in a row. The 
initial weights are assigned to one of the added edges (bold 
edges) that join one of the original nodes and one extra node, 
(d) The original nodes are duplicated preserving their neigh- 
bors. Additionally, the duplicated nodes are interconnected 
by an added edge carrying a zero weight, (e) A MWPM is 
computed: bold edges are matched and dashed ones are un- 
matched, (f) Reconstruction to the initial lattice taking the 
MWPM result into account. The thickly denoted path and 
loop represent the minimally weighted configuration. 



in Fig.[2jc). Therein, one of the two edges connecting an 
additional node to an original node gets the same weight 
as the corresponding edge in the original graph G. The 
remaining two edges get zero weight. 

(3) Subsequently, the original nodes i € V (i.e., not 
the two extra nodes) are "duplicated", i.e. i — > 
along with all their incident edges and the correspond- 
ing weights. For each of these pairs of duplicated nodes, 
one edge {ii, 12} with zero weight is added that connects 
the two nodes i± and i%. The resulting auxiliary graph 
G(V"a, Ea) is illustrated in Fig.J^d) without showing the 
edge weight assignment for reasons of clarity. Note that 
while the original graph is symmetric, the transformed 
graph is not. This is due to the details of the mapping 
procedure and the particular weight assignment we have 
chosen. A more extensive description of the mapping (in 
terms of minimum-energy domain wall calculations for 
the 2D Ising spin glass) can be found in Ref . [IB] ■ 

(4) A MWPM on the auxiliary graph is determined 
via exact combinatorial optimization algorithms |40j . A 
MWPM is a minimum- weighted subset M of Ea, such 
that each node contained in Va is met by precisely one 
edge in M, This is illustrated in Fig.[2je), where the solid 
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edges represent M for the given weight assignment. The 
dashed edges are not matched. Due to construction, the 
auxiliary graph consists of an even number of nodes and 
due to the fact that pairs of duplicated nodes are con- 
nected by additional edges (see step (3)), it is guaranteed 
that a perfect matching exists. 

(5) Finally it is possible to find a relation between 
the matched edges M on Ga and a minimally weighted 
configuration consisting of loops and one path in G by 
tracing back the steps of the transformation. As re- 
gards this, note that each edge contained in M that 
connects an additional node (square) to a duplicated 
node (circle) corresponds to an edge in G that is part 
of a loop or path, see Fig. [2|L). More precisely, there 
are always two such edges in M that correspond to one 
loop/path segment on G. All the edges in M that con- 
nect like nodes (i.e. duplicated-duplicated, or additional- 
additional) carry zero weight and do not contribute to the 
minimally weighted configuration and correspondingly do 
not belong to loops or to the path in the original graph. 
Once such a configuration is found, a depth-first search 
[371 ES] can be used to identify the path and to determine 
its geometrical properties. For the weight assignment il- 
lustrated in Fig. [2|a), the path features u> p — —4 and 
length £ — 3. 

It is important to emphasize that the result of the cal- 
culation is a collection C of loops and one path, such that 
their total weight, and consequently the configurational 
energy E, is minimized. Hence, one obtains a global 
collective optimum. This implies that the resulting 
path is not necessarily the minimally weighted path 
on the lattice: due to the no-crossing condition, loops 
might affect its precise location. While all loops that 
contribute to C possess a negative weight, the weight 
of the path might be positive as well. The set of loops 
might be empty, but the addition of the two extra nodes 
during the first step of the mapping outlined above en- 
sures that the solution of the auxiliary minimum-weight 
perfect matching problem always yields (at least) a path. 

Note that the choice of the weight assignment in step 
(2) is not unique, i.e. there are different ways to choose 
a weight assignment that all result in equivalent sets of 
matched edges in the transformed lattice, corresponding 
to the minimum-weight configuration in the original lat- 
tice. Some of these weight assignments result in a more 
symmetric transformed graph (see, e.g. [37]). However, 
this is only a technical issue that does not affect the re- 
sult. 



III. RESULTS 

As pointed out in the introduction, we aim at finding 
an answer to the question whether paths in the NWP 
model can be described in terms of SLE at the critical 
point, and, if so, which value of k characterizes these 
paths. The critical point p c = 0.3413(9) has been es- 
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FIG. 3: Estimation of the fractal dimension df by means of 
a simple fit to the power law data. The inset illustrates the 
estimation using local effective exponents. Both estimates go 
well together within the error bars. 



timated in the same way as in Ref. [24] with slightly 
higher precision. Principally SLE curves can be calcu- 
lated by Eq. [I] This stochastic differential equation does 
not provide a unique curve for a given diffusion constant 
k, but provides a set of curves that exhibits particular 
properties in a statistical sense. These properties can be 
utilized to approach the question whether a given ensem- 
ble of curves is a potential candidate for showing SLE. 
The first property takes the fractal dimension df into ac- 
count. The fractal dimension can be defined from the 
scaling of the average length (I) of the paths as a func- 
tion of the system size L, according (£) ~ L d > ' . In the 
context of SLE there is a relation between the diffusion 
constant and the fractal dimension [3H [32] : 



K = 8(d f - 1), 



K < 8. 



(4) 



The fractal dimension is fixed at df = 2 for n > 8. A 
simple fit to the power-law data shown in Fig. [3] yields 
the exponent df = 1.283(2). 

An estimation by extrapolation using local effective 
exponents [33] yields df — 1.27(9). Both estimates go 
well together within the error bars. Each data point is 
generated by means of 12800 realizations of the disorder. 
According to Eq. Q, the more precise estimate yields 
k = 2.26(2). 

Similar to [14j [44H47] we examine predictions of 
Schramm's left passage formula (SLPF) in order to ob- 
tain a second estimation for k. As proven in [48J, 
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FIG. 4: (a) Illustration of a minimum-weight configuration consisting of loops (gray) and one path (black) in a lattice of size 
257 x 64. Due to the "Manhattan-structure" of the underlying lattice, the semicircle (R = 32) assumes the shape of a triangle 
(light gray). The boundaries in vertical (bottom-top) direction (dashed) are free and the ones in horizontal direction (solid) 
are periodic. The edge weights are taken from a "Gauss-like" distribution shown in Eq. Q featuring p = 0.3413. The path is 
forced to connect the horizontal boundaries as described in the text. The minimum- weight configuration has to minimize Eq. 
Q. (b) Sketch of a lattice that approximates the upper half plane to check SLPF predictions at some certain places located 
on the semicircle. The checkpoints marked by rhombuses lie on the left and the other ones (crosses) on the right of the path 
(black). For every realization of the disorder, the predictions of SLPF are examined independently for each checkpoint. cj> 
states the angle between the ordinate and considered checkpoints in the upper half plane. 



yields the probability that a curve which can be described 
in terms of SLE will pass to the left of a given point 
z = x+iy in the upper half plane provided that the curve 
links the origin to infinity. 2F1 denotes the Gaussian hy- 
pergeometric function and T the gamma function. P K (z) 
does not depend on the distance between the origin and 
point z. Thus the ratio between Re(z) and Im(z) can 
be replaced by a function of an angle. A proper angle 
4> G [— 7r/2,7r/2] between the imaginary axis and the ref- 
erence point z is located in the origin of the upper half 
plane as depicted in Fig. [4] (b) [H] . Thus the ratio be- 
tween x and y is given byz/y = tan0. A curve 7 will 
be considered as passing to the left of a given point z if 
there exists a continuous path between z and the posi- 
tive real axis that does not intersect 7 [15] . To support 
intuition, Fig. [5] shows SLPF for a particular value of n 
(at this point the reference to the figure is only meant to 
illustrate the general functional form of the SLPF, the 
precise value of the parameters p and k and the results 
implied by them are discussed below). 

Note that Eq. ^ holds in the upper half plane only. 
Subsequently we consider two different lattice geometries 
that aim to approximate the upper half plane. 

On the one hand we approximate the upper half plane 
directly by constructing a rectangular square lattice fea- 
turing appropriate boundary conditions. On the other 
hand we design a square lattice exhibiting a circular con- 
tour and regard it as unit disc E in the complex plane. 
The unit circle can be mapped conformally into the up- 
per half plane by a simple mapping called Cayley function 

g : E — >■ H, z^ii±^. (6) 
1 — z 

The mapping provides g(—l) = and g(l) = 00, thus 
a curve in E linking — 1 + Oi and 1 + Oi is mapped by 
g as required. For a given realization of the disorder, 
the ground state in the unit disc and the corresponding 




FIG. 5: Comparison between the measured probabilities p(4>) 
and the predictions of SLPF P K (cj>). The left passage probabil- 
ity p(<j>) has been measured for various points on a semicircle 
(R — 128). cj> denotes the angle that is illustrated in Fig. [4] 
(right). The computational simulations have been performed 
at the critical point in a lattice of size 1025x256 using 51200 
realizations of the disorder. The chosen diffusion constant 
k* =3.337 provides the finest agreement between Pk,(4>) and 
p(4>) (cf. Fig. |7j. The inset illustrates the deviation of the 
measured left-passage probability from the SLE predictions. 
For reason of clarity, merely a few calculated data points are 
illustrated. As a matter of fact, p(<j>) has been estimated for 
257 different values of <f>. 



ground state in the upper half plane are illustrated in 
Fig. [6] to support intuition. 

First we consider the rectangular lattice shape. In or- 
der to approximate the upper half plane as accurately 
as possible, the terminal points of the paths in the NWP 
model are chosen as described in Section [TTJ The terminal 
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FIG. 6: (a) Illustration of a minimum-weight configuration consisting of loops (gray) and one path (black) in the unit disc 
with a radius of 40. The terminal points of the path are fixed on — I + Oi and 1 + Oi to ensure that they will be mapped as 
required by Cayley's function. All points which Schramm's left passage formula is checked for are located on the vertical solid 
line (deep gray). 

(b) The same minimum-weight configuration as in (a) after the conformal mapping by Cayley's function has been performed, 
so the path connects the origin to infinity in the upper half plane. The dashed lines indicate that the lattice is not displayed 
completely. Due to the mapping, the original, vertical line containing all checkpoints becomes a semicircle. 



point located at the central node of the lower boundary 
is considered as the origin of the coordinate system. Fur- 
thermore, all considered lattices are on a scale of 4 : 1 to 
avoid the influence of the finite size in horizontal direc- 
tion. An example of a minimally weighted configuration 
in such a lattice is shown in Fig. [4] (a) for a lattice of 
size 257 x 64. Although SLPF applies for all points in 
the upper half plane, its predictions are often checked 
for properly chosen points only. Aside from k, Eq. |5| 
depends on angle (j> merely, thus it is sufficient to con- 
sider points located on a semicircle, whose central point 
is located in the origin of the upper half plane. The 
distance between the center of this semicircle and ev- 
ery point located on its boundary should be equal. For 
that reason and due to the "Manhattan structure" of the 
underlying lattice, the semicircle assumes the shape of 
a triangle as illustrated in Fig. [4] The path is able to 
cross the semicircle at particular points only. The check- 
points we consider to examine predictions of SLPF are 
located exactly between these possible crossing points. 
Each checkpoint is examined independently for a given 
realization of the disorder, i.e. for each given path. In so 
doing, for a sufficient large number of samples N s , each 
checkpoint will exhibit a measured left-passage probabil- 
ity p{4>) ± i/p(</>)(l— p{4>))/{N s ~ 1) [51] . A comparison 
between these measured probabilities p(<fi) and the cor- 
responding predictions of SLPF provides an estimate for 
k. The minimum of the cumulative squared deviation 

yields the finest agreement between P K ((p) and p((j>), 
where SC denotes a set containing all possible angles 
of the considered checkpoints. The cumulative squared 
deviation as a function of K is illustrated in Fig. [7] for a 
lattice of size 1025 x 256. 

In so doing fsc( K ) an d its error bar A/ sc (k) have 
been calculated for a set of M = 4000 individual values 
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K 

FIG. 7: The cumulative squared deviation /sc*(«) of the 
measured probabilities p(4>) from the corresponding values 
given by Schramm's left passage formula P K {4>) a s a function 
of the diffusion constant k. The semicircle, which contains 
the checkpoints, features a radius i? = 128 in a lattice of size 
1025x256. The simulations (51200 realizations of the disor- 
der) have been performed at the critical point. The squared 
deviation has been calculated for discrete values of k between 
2.5 and 4.5 with step range 0.001. For reasons of clarity, only 
a few error bars are displayed. The cumulative squared devi- 
ation yields the finest agreement between p(4>) and P K {<fi) at 
k* = 3.337 ±0.077. 



of k £ K , where K = {ki, k,2, ■ ■ ■ , km}- If k* denotes the 
diffusion constant where the squared deviation fsc( K ) is 
minimal, its error bar Ak* will be estimated as follows: 
Defining fi = { K eK \ /sc(«*) + A/ sc (/s*) > fsc^)- 
A/sc(«)}j one has Ak* = max{|K* — n\ | kgO}. This 
results in k* =3.337 ± 0.077. One might wonder whether 
the result is affected by different lattice-size effects, which 
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FIG. 8: (a) The diffusion constant k* depending on the ratio 
between the lateral extensions of the lattice. The horizontal 
extension L x is varied in length and the radius R of the semi- 
circle is chosen to be the half of the system size L = L y . The 
estimates of the diffusion constant have been determined at 
the critical point by means of 51200 realizations of the disor- 
der. For reasons of clarity, the error bars are not illustrated. 
Note that all estimates (except for the shortest horizontal ex- 
tensions until 80) are compatible to each other within the 
error bars. 

(b) The diffusion constant k* depending on the system size L 
that denotes the number of nodes in the vertical direction of 
the lattice. The horizontal extension is equivalent to 4L+1. 
The radius of the semicircle is chosen to be the half of the 
system size in every cases. The estimates of the diffusion con- 
stant have been determined at the critical point by means of 
51200 realizations of the disorder. 



we consider now. 

As mentioned above, we design the underlying lattice 
on scale of 4 : 1 to approximate the upper half plane. It 
is evident from Fig. [8|a) that this ratio is sufficient. We 
estimate the diffusion constant using the squared devi- 
ation technique described above for different aspect ra- 
tios. In doing so we fixed the number L y = L of nodes 
in the vertical direction and varied the horizontal exten- 
sion L x . It becomes apparent that the diffusion constant 
does not vary significantly above a horizontal extension 
that is four times bigger than the vertical extension of 
the lattice, in particular, all values are clearly different 
from the value corresponding to the fractal dimension. 

In order to receive an impression of the influence of fi- 
nite size effects, we varied the system size L (using the 4:1 
aspect ratio) as depicted in Fig. pjb). Bear in mind that 
the system size L corresponds to the number of nodes in 
the vertical lattice direction. It becomes apparent that 
the diffusion constant estimated by using the squared de- 
viation technique is steady sufficiently atop L = 128. A 
further study, which is not illustrated here, reveals that 
k* does not vary significantly for differently chosen radii 
of the semicircle. 

Thus a good approximation to the upper half plane is 
provided by a lattice of size 1025 x 256 in addition to a 
semicircle with radius R — 128. So we estimated k for 
such a lattice by examining SLPF predictions. As men- 
tioned above, the cumulative squared deviation fsc(^) 



of p(<p) and P K (<j>) at p c yields k* = 3.337(77) (cf. Fig. 
W\. A comparison to the estimate obtained from Eq. Q 
(k = 2.26(2)) reveals clearly that both estimates do not 
agree. 

To ensure that this conspicuous distinction does not 
arise from the approximation of the upper half plane by 
designing a rectangular lattice, we also examine paths 
on the unit disc and check the predictions of SLPF after 
mapping the unit disc (and, consequently, the determined 
ground states) into the upper half plane, as depicted in 
Fig. [6] The radius of the disc is discretized by 128 nodes 
and the terminal points of the paths are set as described 
above (cf.[6^a)). Just like before, the predictions of SLPF 
are examined merely for checkpoints located on a semi- 
circle (cf. Fig. gb)) in an analogous manner. Taking 
12800 realizations of the disorder into account, the dif- 
fusion constant k* = 3.34(14) provides the best match to 
the measured left-passage probability. Since this estima- 
tion of k is equal to the estimation obtained before, it 
becomes apparent that the upper half plane has been ap- 
proximated properly by the rectangular lattice and the 
estimate of K obtained by Eq. ^ differs significantly from 
the estimate considering SLPF as an actual fact. Conse- 
quently this means that paths in the NWP model cannot 
be described in terms of SLE. 

Note that we also compared the measured left-passage 
probability p((j>) and the corresponding predictions of 
SLPF directly. From Fig. [5] it is evident that p((f>) and 
-P3.337(</ > ) go well together at almost all value of <f>. Al- 
though the paths of the NWP model do not show SLE 
properties, the predictions of P3. 337(0) appear to apply 
for this model at the critical point. 



IV. CONCLUSIONS 

We examined whether paths in the negative-weight 
percolation model fall into the classification scheme of 
SLE. In order to address this question, we studied their 
geometrical properties and compare them to predictions 
that arises from the SLE theory. First we determine the 
fractal dimension dj — 1.282(2) by measuring the average 
path length as a function of the system size. The relation 
K = 8(df — 1) provides an estimate for the diffusion con- 
stant k — 2.26(2). In order to get a second estimate for 
k, we examined Schramm's left passage formula, which 
states the probability that a curve featuring particular 
boundary conditions and showing SLE passes to the left 
of a given point in the upper half plane. For this purpose 
we designed lattices featuring different shapes in order to 
approximate the upper half plane as sufficiently as possi- 
ble. A comparison between the predictions of Schramm's 
formula and our numerical results provides the estimate 
k* = 3.34(8). Because of the conspicuous distinction be- 
tween both estimates of ft, it became apparent that paths 
in the NWP model cannot be described in terms of SLE. 
We propose two possible reasons for this failure of SLE 
for the NWP. First, the minimum- weight path from the 
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center bottom to the top might be affected by negative 
loops, which arise due to the minimality criterion. The 
edges occupied by these loops are not available to the 
minimum-weight path. Nevertheless, as visible from Fig. 
[I] these additional loops are small and somehow rare, 
hence, the effect might be negligible. Note that algorith- 
mically excluding the loops from the minimization pro- 
cedure will make the problem algorithmically hard, i.e., 
the running time increases exponentially with the system 
size, such that only very small instances can be solved. 

Second, the way the minimum-weight path is selected 
refers to a global minimization criterion, hence it can- 
not grow locally. Loosely speaking, the path can af- 
ford to cross a region containing few negative edges if 
it arrives to a region which exhibits a high concentration 
of negative edges as compensation. This might violate 
the so-called "Markov-property" of SLE processes [18]. 
On the other hand, this argument applies only to single- 
realizations and might be not relevant for the stochastic 
ensemble, which is what the SLE property describes. In 



fact, domain- walls in T = 2d spin glasses are also ob- 
tained by a global optimization procedure, and they seem 
to obey SLE. Thus, the reason for the failure of the SLE 
property for NWP is still not totally clear to us. 
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